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Abstract 

Radial pulsations of neutron stars and strange quark stars with nu- 
clear crust are studied. The avoided crossing phenomenon occurring for 
the radial modes is found and discussed. Neutron star models are con- 
structed using a realistic equation of state of dense matter and strange 
star models using a phenomenological bag model of quark matter. The 
eigenfrequencies of the three lowest modes of linear, adiabatic pulsations 
are calculated, using the relativistic equations for the radial oscillations. 

1 Introduction 

Radial oscillations of dense stars (neutron and strange stars) were studied by 
many authors (see eg. Glass & Lindblom ( jl983| ), Vath & Chanmugam ( p92j ) 
and references therein). 

In the present paper we study the radial pulsations of neutron stars and 
strange quark stars with nuclear crust in the framework of general relativity. 
Our aim is to present an interesting phenomenon - avoided crossing of this ra- 
dial modes. Avoided crossing manifest itself as a rapid change of the frequency 
of the certain oscillatory mode when stellar configuration gradually changes. 
At this point the frequencies of two consecutive modes are very close and star 
"switches" the oscillatory properties from one to another. This effect was pre- 
viously noticed in the analysis of the nonradial oscillations in "ordinary" stars 
(Aizenman, Smeyers and Weigert ( [ 1977 1 ) , Christensen-Dalsgaard ( ||1980fl )) and 



has been recently discussed by Buchler, Yecko and Kollath ( [ 1997 1 ) in the case 
of radial pulsations of the classical variable stars. For neutron stars the avoided 
crossing phenomenon was considered by Lee and Strohmayer (| 1996| ) for non- 
radial oscillations of rotating neutron stars. 

The avoided crossings of radial modes in neutron stars was mentioned in our 



recent paper (Gondek, Haensel and Zdunik ( [ 1997 1 , hereafter Paper I) also in 
the case of hot protoneutron star. 



In this paper we study the avoided crossing phenomenon for two models of 



dense matter: the low temperature limit of the Lattimer and Swesty ([1991 



equation of state of the neutron matter and the strange quark matter described 
by the simple bag model. 

The plan of the paper is as follows. In Section 2 we describe equations of 
state for which we found the avoided crossing phenomenon. We calculate here 
also the adiabatic index of the matter which is crucial for the existence of the 
avoided crossings. In Section 3 we present the formulation of the problem of 
linear, adiabatic, radial pulsations of stars. We discuss here also the method of 
determination separately the oscillatory properties of the inner part of a star 
and the outer regions (envelope). Numerical results for the eigenfrcqucncies of 
the lowest modes of radial pulsations of stars are presented in Section 4. Finally, 
Section 5 contains a discussion of our results and conclusions. 



2 Equation of state and adiabatic indices 



The starting point for the construction of our equation of state (EOS) for the 
neutron star models was the model of hot dense matter of Lattimer and Swesty 
( 1991 ), hereafter referred to as LS. Actually, we used one specific LS model, 



corresponding to the incompressibility modulus at the saturation density of 
symmetric nuclear matter K = 220 MeV. We use LS model of matter in the low 
temperature limit, which can be treated as a zero-temperature T = model. 

The detailed discussion of the EOS and equilibrium conditions is presented 
in Paper I. 

As a second example we consider the strange stars with nuclear crust. The 
interior of this star is build of strange quark matter, containing nearly equal 
number of u, d and s quarks. This matter is described in the framework of the 
bag model with the v alue of bag constant B = 60 MeV&n. 3 (Wittcn ([1984]), 
Farhi & Jaffc ( [jl984| ), Haensel, Zdunik & Schaeffer fltt986| )). The crust of a 
quark star, at densities below neutron drip point (/?nd = 4.3 • 10 11 gem" 3 ), is 



represented by the BPS equation of state (Baym, Pethick & Sutherland | 1972| ). 
In Fig. [l] we show our equations of state for nuclear and quark matter. 



2.1 Adiabatic indices 

The adiabatic index within the star, T, plays a crucial role for oscillatory prop- 
erties of the star. 

The oscillation frequency depends sensitively on the value of the adiabatic 
index and on the shape of the T(p) function. The increase of L means the 
growing stiffness of the matter leading to the effect of the "wall", setting the 
bounds for the regions of oscillatory motion. 

To calculate properly the oscillations of the star one has to know the rel- 
evant index, governing linear perturbation of the pressure due to the density 
perturbation. This index will be denoted by r osc . 
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Figure 1: Pressure versus baryon density for our models of dense matter. The 
solid curve corresponds to cold catalyzed nuclear matter. The dotted curve 
presents the EOS for strange quark star with nuclear crust. 
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Figure 2: Parameter r osc versus pressure, for the nucleon matter. 

In the case of a sufficiently cold matter, the reactions between matter con- 
stituents are so slow, that the matter composition remains fixed (frozen) during 
perturbations in the dynamical timescale, because T roa ct 3> 7d yn - in such a case, 
the appropriate adiabatic index is rf rozon i.e. 



where s, Y e correspond to the equilibrium model. Fixed value of Y e during oscil- 
latory motion results in freezing of the weak interaction processes (/3- reactions) . 

The detailed discussion of the relevant adiabatic indices T osc at different 
physical conditions (hot protoneutron star with or without trapped neutrinos) 
was presented in Paper I. 

The dependence of T's on the pressure in the stellar interior, for our models 




(1) 




Figure 3: Parameter r osc versus pressure, for the strange quark matter. The 
left figure in logarithmic scale visualizes the huge jump in T at the crust-strange 
core boundary. In the right figure we presented T in linear scale in the region 
of core oscillations. For comparison we show T for nuclear matter (dashed line) 



is displayed in Figs, [|. In this figures (and also some next figures) we choose 
pressure as an independent variable because pressure is continuous through the 
star and describes very well some stellar parameters (e.g. mass above given 
radius). This is not the case of density p, which exhibits large density jump at 
the crust-core boundary of quark stars. 



3 Linear adiabatic radial pulsations of stars 

The stellar configurations in hydrostatic equilibrium are determined by the into 



gration of Tolman-Op pcnh cimcr-Volkoff (TOV) equations (Tolman [1939], Op- 
penheimer & Volkoff [|l939| ). 



The mass-radius relation for the neutron and quark star models is shown in 
Fig. |U In Fig. [| we show mass vs. central pressure for our stellar configurations. 
The points on the M(R) curve for neutron stars correspond to the stellar models 
for which oscillatory properties were studied in details and presented in Fig. |^ 
and discussed in Sec. 4 . 

To find eigenfrequencies u> we solve the equations governing infinitesimal 
radial adiabatic stellar pulsations in general relativity derived by Chandrasekhar 
(| 1964| ), and rewritten by Chanmugan ( 1 1977 ] ) in a form, which turns out to 



be particularly suitable for numerical applications. Two important quantities, 
describing pulsations, are: the relative radial displacement, £ — Ar/r, where 
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Figure 4: The gravitational mass versus stellar radius for static models of neu- 
tron stars (solid curve) and quark stars (dotted curve). The points on the M(R) 
curve for neutron stars correspond to the stellar models for which oscillatory 
properties were studied in details and presented in Fig. pi 
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Figure 5: The gravitational mass versus central pressure for static models of 
neutron stars (solid curve) and quark stars (dotted curve). The points on the 
M(R) curve for neutron stars correspond to the stellar models for which oscil- 
latory properties were studied in details and presented in Fig. pi 
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Ar is the radial displacement of a matter element, and AP - the corresponding 
Lagrangian perturbation of the pressure. These two quantities are determined 
from a system of two ordinary differential equations which can be reduced to the 
one second order (in £), linear radial wave equation, of the Sturm-Liouville type 
with lo 2 as the eigenvalue of the Sturm-Liouville problem. As a result for a given 
stellar configuration, we get a set of the eigenvalues lu 2 < lo\ < ... < uj 2 < .... 
with corresponding eigenfunctions £o, £1, £m where the eigenfunction £ n 



has n nodes within the star, < r < R (see, e.g., Cox [198C]). The detailed 
form of the equations and solving method is presented in Paper I. 

To solve linear radial wave equation one needs two boundary conditions. 
The condition of regularity at r = requires 

(AP) center = -3(^rp) center . (2) 

The surface of the star is determined by the condition that for r — > R, one has 
P — > 0. This implies 

(AP) surface = (3) 
To study the oscillations of the core end envelope separately we divide the 
star into two parts: core from the center down to the boundary density pb and 
envelope with density smaller than pb- In the case of neutron matter the density 
Ph is defined by the steep change in T (Fig. ||) and is equal to ~ 10 14 gem -3 . 
The corresponding pressure is equal to Ph — 10 32 dyncm -2 . For quark stars 
the crust-core boundary is very well defined and connected with huge density 
jump at the pressure Pb = 10 30 dyncm -2 . At this boundary the density of the 
crust is equal to Pnd = 4.3 • 10 11 gem -3 and the density of quark core is slightly 
above 4P/c 2 = 4.28 • 10 14 gem" 3 . 

We have to apply the boundary conditions appropriate for core or envelope 
pulsations. In the investigation of the core oscillations we stop the integrations 
of the linear radial wave equation at the Pb with the boundary conditions cor- 
responding to the "free" envelope oscillation with the frequency defined by the 
central core - we treat the envelope as a mass M cnv laying on the central core 
and simply moving in a way governed by this central region. This assumptions 
leads to the boundary condition in the form: 



(4) 



AP / (GM lo 2 R 3 \{ 2GM\ l V 
— = -{ 4+ {rc 2 - + Gull 1 - -R^) 
The oscillation of the envelope we can define by the boundary condition: 

£0%) = (5) 
which means that central core do not move (in other words £(p > pb) = 0). 



4 Eigenfrequencies 

As a result of the numerical integration of linear wave equation with bound- 
ary conditions appropriate for three situations (Eqs. (|[ ||, [| ||)) we obtain 
eigenfrequencies of oscillations for each stellar configuration. 
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Figure 6: The eigenfrequencies of the three lowest modes, versus stellar mass for 
neutron star's radial oscillations. The solid curve corresponds to the oscillations 
of the star as a whole, the dotted curve describes oscillations of the crust of a star 
and dashed curve - oscillations of the core. The points on the M(R) curve for 
neutron stars correspond to the stellar models for which oscillatory properties 
were studied in details and presented in Fig. || 
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Figure 7: The cigcnfrcquencics of three lowest modes, versus stellar mass for 
quark star's radial oscillations. The solid curve corresponds to the oscillations 
of the star as a whole, the dotted curve describes oscillations of the crust of a 
star and dashed curve - oscillations of the core. 
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In Figs. ||, [?| we plotted the eigenfrequencies of the three lowest n = 0, 1, 2 
radial modes versus stellar mass for neutron and quark star oscillations respec- 
tively. The solid line corresponds to a whole star pulsations with boundary 
conditions at the center and stellar surface (Eqs. |3|)). The avoided cross- 
ing phenomenon is clearly visible and manifest itself as abrupt, nearly stepwise 
changes in frequencies of the consecutive (neighbouring) oscillation modes. This 
effect occurs also in the case of hot, isothermal models of protoneutron stars 
(Paper I). 

This stepwise changes of oj n are due to the change of the character of the 
standing-wave solution for the eigenproblem. Namely, at the avoided crossing 
point the solution changes from the standing wave localized mainly in the outer 
layer of the star, to that localized predominantly in the central core. This is 
very clearly visible from the comparison with solutions of oscillatory equation 
for core and crust pulsations separately. The solid line of the frequency of 
stellar pulsations nearly coincides with either "core" or "crust" solutions. With 
increasing mass of the star the subsequent modes u> n (M) choose between the 
crust or core solution to fulfill the condition uj^ < u\ < ... < w„ < Because 
of the different dependence uj n (M) for crust and core oscillatory modes this 
leads to the avoided crossing effect. 

The crust region in the stars under consideration is relatively small. Thus 
the oscillations of the crust could be very well described neglecting the mass of 
the crust in the gravitational force. In the Newtonian limit and for constant 
T = r crust throughout the stellar crust the wave equation for crust oscillations 
can be reduced to the form of Bessel equation of the order of v defined by the 
formula: 

r = i + - (6) 

v 

The resulting oscillatory frequency should then be proportional to the gravity 
at the stellar surface 

1 3u,n 




2,7(l + I)Vft y (7) 

where j Vtn denotes the n-th zero of the Bessel function J„ (z) of the order of v and 
g = Sgf is gravitational acceleration at the stellar surface. This approximation 
agrees very well with our exact numerical results for neutron and quark stars. 
Especially for strange quark stars where T is nearly constant and equal to 4/3 
in the large part of the crust (Fig. ^) the difference between true, numerical 
results and approximation by the formula (Q) with v = 3 is smaller than 5%. 
In the case of neutron stars due to the changes of T through the outer layers 
the identification of subsequent oscillatory modes with Bessel functions is less 
accurate, but the linear dependence of uj n on g in formula (]?]) holds very well. 

To confirm our conclusions we have studied in details properties of the three 
lowest radial modes for some stellar configurations. The parameters of this con- 
figurations and eigenfrequencies w n are presented in Figs. by dots A-F. 
In Fig. U we plotted the energy density of oscillations calculated per pressure 
gradient in the star i.e. e osc (P) dP is equal to the energy due to the oscillatory 
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Figure 8: The oscillatory energy (e osc ) vs relative pressure (P(r)/P cen tr) f° r 
three lowest radial modes, n = 0, 1, 2 for neutron star. The subsequent rows 
(marked by letters A-F) correspond to stellar models with growing central den- 
sity (from A to F) marked by dots in Figs. ||, ||, q. The energy is normalized to 
unity i.e. J e osc d(P/P CC ntr) = 1< The function P(r)/P cen tr approximates very 
well the relative mass from the stellar surface (M — m{r))/M . 
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Table 1: Identification of radial oscillatory modes of neutron star (see Fig. 



propagation 


fundamental 


first 


second 


region 


mode 


overtone 


overtone 


crust 


AO BO CI E3 


Al B2 




core 


Bl CO DO EO FO 


D2 El Fl 


F2 


mixed 


A2 C2 Dl 





motion of the matter in the sphere r(P) — r(P + dP). This function is pro- 
portional to £n( r )- As an independent variable we choose in Fig. || P(r)/P centT . 
For the outer layers of the star this function is proportional to the relative mass 
above the radius r i.e.: 

P(r)/P ccntI ~ (M - m(r))/M (8) 

with proportionality constant of the order of unity. The integral of the energy 
in Fig. [|is normalized to unity i.e. J e osc d(P/P C entr) = 1. 

The Fig. || very well visualizes where the mode is localized. We can also 
identify modes by the shape of £ or energy of oscillation. The three lowest 
radial "core modes" are characterized by the functions presented in row F in 
Fig. ||. The stellar configuration F is the most massive and crust plays very little 
role in oscillatory properties (see also Fig. ^|). The two lowest "crust modes" 
correspond to the columns n — and n = 1 in row A of Fig. g[ For the 
low massive model A the crust is relatively large and determines frequencies 
for fundamental mode and first overtone. Second overtone n = 2 is strongly 
contaminated by the core oscillations (Figs. || and ||). Fig. || allows for very 
clear division of oscillatory modes by the real pulsational properties and not the 
number of nodes. This is presented in Table 1, where for example Al denotes 
configuration A, number of nodes n = I. 

The identification of the radial modes in Table 1 is consistent with the Fig. |[ 
The specific oscillatory mode for stellar configurations A-F corresponds to the 
crossection points of the vertical lines at points A-F and the curves ui n (M) 
for given stellar configuration. The location of these points very well define 
type of oscillatory motion. The nearly constant parts of the w a (M) functions 
correspond to the core oscillations while regions of increasing u> n (M) describes 
crust oscillations. The region close to avoided crossing points is characteristic 
to mixed modes: star pulsates in crust as well as in the core. 

5 Discussion and conclusions 

The avoided crossing phenomenon in the radial pulsations of dense stellar objects 
was presented. Both strange quark stars with crust and neutron stars could be 
subject to this effect. 
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The main reason for the abrupt changes of the oscillatory properties of the 
star as one change a little stellar configuration is the very well defined division 
of the star into two regions: inner core and outer crust. This two parts have dif- 
ferent equation of state and oscillatory properties. The latter have been studied 
in this paper by the separation of the outer and inner regions by appropriate 
boundary conditions. Resulting oscillatory frequencies of crust and core pulsa- 
tions depend on the stellar configurations (e.g. mass of the star) in completely 
different manner. The functions w n (M) crust are rather steeply increasing while 
w n (M) core are nearly constant (neutron stars) or decreasing (quark stars). 

The oscillatory properties of the star could be characteristic either to the 
crust or core pulsations. The frequency spectrum depends on the mass of the 
star. The core pulsations determine the oscillatory properties of the lower order 
radial modes for relatively massive neutron or strange stars. The star could 
be subject to the crust pulsations for sufficiently high frequency or low stellar 
mass. For given mass the fundamental mode of crust oscillations uo(M) crvLE t 
define the frequency below which the star pulsates in the core (Figs. ^, 0). 

The avoided crossing phenomenon is strongly related to the changes of com- 
pressibility of the matter throughout the star, which is described by the shape 
of the r osc (yo) function. The increase of the stiffness of the matter outward leads 
to the maximum of T osc (p) close to the boundary between core and crust. In 
the case of strange quark matter this is the consequence of the the fact that it 
is self-bound at the very high density and r — > oo as P — > 0. In neutron matter 
described by LS equation of state the symmetry energy leads to the increase of 
Y e inward and in this case more symmetric matter is softer. The same effect 
could be seen for the field model of dense matter considered by Diaz Alonso 



and Ibahez ( 1 1985 ). The function T(p) for their EOSs (Fig. 9 of their paper) 



resembles LS result at low temperature limit. 
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